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Abstract 



Ph. 



For the one-dimensional Schrodinger equation, some real intervals with no eigenvalues (the spectral gaps) 
may be obtained rather systematically with a method proposed by H. Giacomini and A. Mouchet in 2007. 
The present article provides some alternative formulation of this method, suggests some possible generaliza- 
. tions and extensively discusses the higher-dimensional case. 

> 

H. Giacomini et A. Mouchet out propose en 2007 une methode permettant d'obtenir des trous spectraux d'un 
operateur de Schrodinger a une dimension, c'est-a-dire des intervalles ne contenant aucune valeur propre. Le 
present article introduit une formulation differente de cette methode, suggere des generalisations possibles 
et traite de fagon exhaustive le cas de plusieurs dimensions. 

> 

in Introduction 

I A spectral gap (or eigengap) of a self-adjoint operator is a closed real interval to which no eigenvalues belongs. 
In [2 we have presented a systematic method for finding gaps in the discrete part of the spectrum of a one- 
dimensional non- magnetic Schrodinger equation with a potential V{x). When V is a polynomial half-bounded 
from below, the boundaries of the gaps are given by the real zeroes of a family of polynomials whose degree D 
may be arbitrary large. The construction of these polynomials is provided by an explicit and straightforward 
algorithm. For still not understood reasons, it happens that in every case we have considered, our method works 
surprisingly well when compared to numerical computations: when increasing D the more and more numerous 
^ ' intervals we compute resolve the spectrum from below (i.e. the lowest eigenvalues are separated by at least one 
I gap) and the infima of each interval seem to converge quickly to the eigenvalues. For the moment, we have no 
clue to understand these two phenomena and our method comes without any estimation of the distance between 
the gaps and the spectrum. 

Being as local as possible (no computation of integrals is required), our method differs strongly in spirit 
from other spectral approximations like the Rayleigh-Ritz variational methods or the Rayleigh-Schrodinger 
perturbation methods. Since generically the spectrum cannot be determined exactly, finding such gaps may 
offer a valuable piece of spectral information, complementary to the information obtained by other methods. 

One natural issue is to try to extend our method to multi-dimensional systems. Then, the nature of the 
spectrum depends on the integrability properties of the system (among the vast literature on this subject see for 
instance [3]). Qualitatively the statistical distribution of the eigenvalues exhibits different correlations according 
to the nature of its symmetries. For instance, for non-integrable systems, the two point correlation function of 
the discrete spectrum exhibits a so-called "level repulsion" because, unlike what occurs in integrable cases, the 
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probability of finding two successive eigenvalues whose distance is s vanishes when s tends to zero. Therefore we 
expect that this dichotomy between integrable and non-integrable cases should somehow appear in any general 
method for finding spectral gaps. However, most unfortunately, we have not been able to generalise our strategy 
in higher dimensions. Although there was a priori no obstacle in sight to such an attempt, it happened that 
the origin of the obstruction came from very subtle arguments that are deeply hidden. One aim of this note is 
to explain (in section [2]) this negative result with the hope that it may help to find out some way to bypass the 
pitfalls or, at least, to help avoiding the same tracks. 

Our second aim is more optimistic but still rely on speculative grounds. After recalling the main ingredients 
of our method in section 11.11 I will introduce a systematic algebraic approach which at first sight seems to 
rephrase in a more elaborate way what we have done in [2]. However by associating with the Schrodinger 
equation the closed algebra of differential operators that will be introduced in §§ II. 2[ II. 3[ we can easily guess 
a fruitful strategy to deal with spectral problems associated with more general (Id) linear equations — for 
instance of order larger than two — or even with non-linear equation like the Gross-Pitaevskii equation (g is a 
real coupling constant) 

ifi" ^ (V - E)ip + gip\ (1) 



1 The one dimensional case: an algebraic approach 
1.1 The principle of the method 

The stationary one-dimensional non-magnetic Schrodinger equation can be written as follows 



dV 

da; 



■j^ = 2{V{x) - E)^ . (2) 



where we will take the real potential V to be smoothly defined on R. The real E will be an eigenvalue 
whenever the real function Lp is square integrable on R. The key idea of our method is to construct, for a given 
integer iV ^ 1, a real function JN{ip' , (p, x, E) (the prime stands for the derivative A/ Ax) such that 

d 



(i) —(Jn{v\x),^{x),x,E)) ={^{x)f Fn{x,E) 



(ii) lim j(^ip'{x),ip{x),x,E) — 0. 

The real function of both the spatial coordinate x and the energy E is (f -independent and is obtained from 
the potential V and its derivatives. For instance we will find the following expressions: 

Fi=-a'^ + 2{V -E)a^, (3a) 

F^^]^4' -A{V-E)a'^-2V'a^, (3b) 

F,=-\ + ^{V-E)a'^ + f V a', + {V"-&{V- Ef) a, , (3c) 
+ (^-\v"'+^^V'{V~E)^ a,. 



(3d) 



where oq is any smooth function such that 



lim ao{x)\ipf = (4) 

\x\^oo 
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which is not very restrictive since we know using semi-classical arguments that Lp itself is exponentially decreasing 
at infinity [T] chap.fO]: 

^{x) ^ e-/'^'V2(n-)-^^)d- (5) 

Condition (i) is the cornerstone of our method and, before justifying how it can be obtained (we will see 
that condition (ii) is not so restrictive), let us first explain how gaps in the spectrum may be obtained. 

When the conditions (i) and (ii) are simultaneously fulfilled, an immediate consequence is that the inte- 
gral j^"^ {,^{^)) F]y{x,E)dx vanishes. This implies that, if E is truly an eigenenergy, the function x 

(^(p{x)^^ Fff{x, E) should change its sign. If N is even, we obtain a 1^9- independent condition: for any fixed 
energy x 1— > Fm{x,E) must change its sign on the real axis. For such a one-dimensional problem, and for a 
given N , we still can choose in a wide continuous set of smooth functions on the real axis because we have a 
lot of freedom in choosing ap. A forbidden value oi E {i.e. E cannot be an eigenenergy) is obtained if we are able 
to chose ao such that remains positive on the whole x-axis. Once this property is achieved, it remains stable 
under small perturbations within the set of ag's, for instance by varying the control parameters A on which oq 
may depend, and we obtain a whole interval where no eigenenergy can exist. More precisely, if we introduce 
explicitly the A-dependence in Fjv, the boundaries of the gaps will necessary be given by some solutions of the 
system of equations 

Fn{x,E,\) :^{)- (6a) 
d^FN{x,E,\)=Q- (6b) 
dxFM{x,E,\)^Q . (6c) 



Using the implicit function theorem where 



\dEFN\\dl,FN\y^O, (6d) 



the first two equations define implicitly a;(A) and E{X); these are the conditions for a bifurcation in the zeroes 
of a; I-)- Fj^{x, E, A) to occur. On one side of the bifurcation x >-)• F]sr{x, E, A) has locally a constant sign (and 
therefore the corresponding value E{X) is forbidden) whereas on the other side x F]\[{x, E, A) locally changes 
its sign and E{X) cannot be ruled out from the spectrum. To put it differently, in the {x, E, A), the set of zeroes 
of Fat becomes tangent to the x-space. Then for each value of A, E{X) is a candidate for being the boundary of 
a spectral gap. If this is the case, we can reach an extremal value provided = dE/dX = —d\F]\[/dEFN that 
leads to the equation (|6c|. 

See 2] for an effective implementation of this method and for applications. In the following, we will remain 
at a more formal level and let us start by defining some notations. 



1.2 Algebraic construction of condition (i) and classification of the possible F/v's 

Denote by Vn, the vector space of smooth applications from R to R^+^. Any element a of Vn niay be 
represented by a vector field x 1— >■ {<in{x)) ^^^^ = (ao(x), . . . , aAr(x)) and can be associated in a one-to-one 
correspondence with the homogeneous polynomial of degree N in the two variables 3> and ^I^: 

N 

Pa{^, $, X) «n(2^)*'^""<f" ; (7) 

71=0 

then, it may be used to construct the real function on R defined by 

N 

Pa{v'{x),^{x),x) 1^^^^a„(a;)(^'(x))^-"(^(a:))" . (8) 

ri=0 
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For computations we will distinguish the "total" derivative D of a function P(^ip'{x),(p{x),x) from its partial 



derivatives d^pi , and 



DP{ip'{x),ip{x),x) = v{x)(p{x)d^'P(ip'{x),ip{x),x) +ip'{x)d^P(ip'{x),ip{x),x) +dxP{ip'{x),(p{x),x) (9) 

where the substitution ip" — vip has been made since we suppose that tp fulfills For simplicity we have left 
implicit the i?-dependence in 

v{x)=^2{V{x)-E). (10) 

From its very definition, it is obvious that the set Vn oi homogeneous polynomials of degree N in and ^p{x) 
is stable under D and, moreover, D is represented by a linear operator D in Vn'- 

/ ao \ ( 1 

Nv 2 '■■ 

Q {N-l)v d.j, 3 



D 



ai 

0.2 



a.N-2 



0\ / ao \ 

ai 

02 



{N-n + l)v 71+1 



V 



3v 






N -1 
2v N 

V d,/ 



aN-2 
\ aN J 



(11) 



For each n E {0, • • • ,N}, denote by Q„ the subspace of Vn defined by a„ = and Q„ its complementary 
defined by the direct sum decomposition Vn = Qn ffl Qn- Looking for all the Jn that fulfill condition (i) can 
therefore be interpreted as the determination in Vn of the preimage D~^Qn- In |2j we have shown how to 
straightforwardly compute Jn and obtain Fn but let us propose a strategy based on a more algebraic formalism 
that may be useful as a warming up for higher dimensions. 

From any a, we can systematically reduce the degree in (p' of Pa if we work up to a total derivative. Indeed, 
for any monomial caracterised by n G {0, • • • , N — 1}, we use the identity 



n+1 



n+l 



1 / lN~n-l n+l N n I iN-n-2 n+2 , n 

■a„(p — ; — va,^ip ip ^ + D 



1 



(12) 



a„ ] , (13) 



n + 1 " n + l ' ' ' ' \n+l 

where, again, we have substituted ip" by vtp. The operation that transforms a„(^'^^"(/j" to the two first terms 
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in (|13l) may be linearly represented in "P^r by the reduction operator defined by 



/ ao \ 




/ 

















o\ 


/ ao \ 


ai 






















ai 


a2 




(1 


-N)v 


--2d. 












a2 




dcf 




















a„ 











n-l-N^ 
n — 1 


n 









a„ 

















]_g 








ajv-i 


\ ajv / 
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V 


\ ajv / 



(14) 



For = 1, we have R = ( ? ) • By construction we have a — Ra = for any vector a = (0, . . . , 0, a^r) e Qn', 
by (jl3p . we have a — Ra € ImD for any vector a £ Q„ with n G {0, • • • , — 1} ; therefore, by linearity, we have 
a — Ra G ImD for any a G Vn ■ Moreover, since R is a lower triangular matrix with all diagonal terms but one 

being zero, all the components of R^'^a vanish but the last one Ajv[a] '= (R^a)jv which is a differential operator 
on a of order iV. For instance we have 

Alia] 

A2[a] 

^3 [a] — — -u-o -T g'^"'0 g I' "0 ^ g"i — 2 ""^ ^ 3"^ ^ "'^ ' 



1 



1 



^0 

1 

6«o 



wa. 



^4 [a] =— ao 



1 

24^ 



— ai 
2 ^ 

7 , 

2 „ 

3«ao- 



3«ao 
3 , , 



1 . 



1 

—va 
2 



1 , 

— On 

3 ^ 



1 , 



12 



3-a, 



1 , 

4«3 



a4 



Now if we use the decomposition 



1 = R 



N 



N-1 

(i-R)ER" 

n=0 



(15a) 
(15b) 

(15c) 

(15d) 

(16) 



it can be seen immediately that any a can be uniquely decomposed in R^a G Qn plus a vector in ImD. Translat- 
ing this decomposition into the language of functions and taking for a the vector associated to DKn(^ip' {x), (p{x), x) 
where is any homogeneous polynomial of degree N in (^(p'{x), 93(2;)) , we have shown that there always exists 
a homogeneous polynomial KN(^ip'{x),(p{x),x^ of degree in (_ip' (x), (p{x)^ , such that 



DKN{^'ix),^{x)) = FNix){^ix))^ - DKN{^'ix),^{x)) 



(17) 



and therefore, in order to recover (i), it is is sufficient to choose Jn = Kjs; + Kjsi. The function Fn{x) is 
independent on ip and ip' and is just given by the action of the linear operator A^r on the coefficients of DK^^. 

A priori, we can start with any set of trial functions (a„)„g{i....^jv} to build up our Kj^, then compute Fn 
by computing the A^**^ power of R. Before we try to control the sign of Fn for even A^, the only restriction 
so far on the a's is to preserve (ii): a„ should not increase faster than (p' 



fN- 



at I a; I — >■ 00. However we will 
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now show that our freedom is in fact restricted to the choice of one test function only. In other words, many 
different choices of (a„)„ will lead to the same Fn and therefore will not help to gain any piece of information 
(in particular those leading to an identically vanishing Fn). To put it very qualitatively, Qat is a very thin 
subspace in Vn (of co-dimension N if seen as a vector space on smooth real functions) and the kernel of D is too 
small (given by the solutions of a linear ordinary differential equation of order N) for D~^Qm to decrease its 
codimension. To understand that, let us introduce the projector II„ on the n**^ component of a. Our previous 
construction of condition (i) can therefore be re-written 

L)a = R^ba . (18) 

where a is the element of Vn associated with the function J^. In section [1.31 we will show directly for small N 
that 

R^b{l-tio) = Q (19) 

which has the following consequence: adding to a any vector h = (6n)n whose 6o = will not affect the 
left hand side of (|18l) from which Fm is computed. Therefore Fm depends only on one function, namely a^. 
All the others can be canceled without loss of generality. Actually, if we start with a = (ao,0, ... ,0), we have 

a Da = (oq, Nvao, 0, . . . , 0) — (oq, ai, 0, . . . , 0) and substituting d with a in (jl5p leads straightforwardly 

to ([2]). If we start with a — (0, oi, 0, . . . , 0), we have d =^ Da — (ai, a'^, (N — l)wai, 0, . . . , 0) — (ao, ai, 0, . . . , 0) 
and (fTO]) can be (tediously) checked in the special cases iV = 2, 3 and 4 : by substituting a with a in (jisp . Aat 
identically vanishes. The same remains true for a = (0, 0, 02, 0, ... , 0), a = (0, 0, 0, 03, 0, ... , 0), etc. The next 
section provides a systematic way of proving the last results and ? 13.11 gives a more general argument. 



1.3 The closed algebra 

Let us prove by introducing an algebra of operators that may be represented by their action on or on 
the functions in V^. The operator ^= ^pd^' when multiplied by v allow to formalise the substitution of Lp" 



by vp. We also define = ip'd^. The arrows recall that these operators raise and lower the component of a: 



/ ao \ 

ai 

ajv-2 





Nao 
(iV-l)ai 

3aAr-3 
2aAr-2 



ai 

0.2 

O.N-2 
O.N-1 



( 



ai 
2a2 
3a3 



(TV - l)ajv-i 
Nan 



\ On 


/ V 


On 


-1 / 




\ aN J 


V 





If we introduce the diagonal matrices A-^ 


dcf ,. 

— dia 


g(o,i,. 


..,N), A^ 


= diag(iV,. 


lowering operators defined by the following action on 


any a 








/ ao \ 
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/ ai \ 




ai 
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ai 
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S+ 










O'N-2 




ajv-3 






aN-2 




ajy-l 




ajv-i 




O.N-2 










ajv 




\ aN ) 




\aN-i/ 






\ ajv / 




V / 



(20) 



,1,0) and the raising and 



(21) 



We have = S Ar^, — S'+A^. Within this formalism, the operator ((TT]) is 

b = S^ + d^+vS^ ^ 5+Ax^ + 4 + vS-A^^ , 



(22) 
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with dx '= d\a.g{dxT ■ ■ ^dx)- To implement the commutation rules between and any diagonal matrix 
A = diag(Ao, Ai,...,AAr), we define the diagonal matrices 

A(i) diag(AiA2, . . . , Aw, Ao) ; A^-^) diag(AAr, AqAi, . . . , \n^i) , (23) 

and for any strictly positive integer k, we define A^^) =^ (• • • (A(i))(i) • • • )^^\ A(-'=) =^ (• • • (A(-i))(-i) • • • 

where the superscript occurs k times in the right hand side. We also take A^"^ A. Therefore the -n}^ entry of 
A^^) for any integer k is (AC^))^^^ = A„+fc modulo n+i- Then we have 

S^A^A^^^^S^. (24) 
In particular if we set II_i = IlAf+i = 0, for any n E {0, . . . , N}, 

S^Iln = fl„^i^± . (25) 
We can also express the reduction operator R in terms of the shift operators S. From (jl4p . we have 

R = flN + A + B = flN + C. (26) 

with 

A ~S-Aadx; (27a) 

B -{S-yhv; (27b) 



and 



A„ diagfl,...,^,...,-!- ) ; (28a) 



n+1' N + 1 

h = diagf jV-l...., ^~"~\ ...,^— ,0,— ^) . (28b) 



From we have S Hat — and therefore CHat = 0. Besides, has a unique non-zero element, 
namely {C'^)n,q and then XlArC*^ = (7^. Therefore, the expansion of i?^ = (IItv + C*)^ leads to 

^ Un{1 + C+--- + C^) ; (29) 

= IItv XI ••• . (30) 

words w — (l,m) 
0^\l\ + \m\^N 

The last sum involves all the distinct words w that can be made with the two "letters" A and B whose length 
is \l\ + 1 771 1 between and N (the zero- length word is the identity). Each word w is uniquely defined by two 
multiple indices, i.e. by two sequences I = (^i, Z2, • ■ • ) and m — (jn\, 7712, . . . ) of positive integers that all vanish 
above a finite rank. We define |^| l\ + I2 + ■ ■ ■ and |777,| 7771 -I- 7772 + • • • • Eventually, the identity ([T^ holds 
if we can cancel, for each 77, G {1, . . . A^}, the product 

^^M„ = nw( X i'iB"ii'^B"^---)(4 + w5-A^-H5'+A^)ll„ . (31) 



words w- 



0^\l\ + \m\^N 

To compute this sum, we start by rearranging each word in order to shift, with the help of ([24|), all the S'^'s 
involved in A and B to the right: 

^;i^mi^/2^m2 . . . ^ A^(S'~)'''+^'"' (32) 
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with the diagonal matrix 

" « ' « ' 

li factors mi factors 

^(-/i-2mi-l) ^ ^ ^ ^(-;i-2mi-/2) Qh ^(-il -2mi -(2 -2) ^ ^ ^ ^(-ii -2mi -(2 -2m2 ) ^ma , , , j-gg-j 
^ a ^ ^^_b ^ 

/2 factors 7712 factors 

Now each word contributes to ([?T|) via three terms 

A„nAr(^-)'"+"'"'~'n„^-5+A^ = A^(5-)'"+"'"'"'n„<5Ar_|z|.2|m|+i.„A^ ; (34a) 

A„nAr(^-)'"+"'"'n„4 = A^(5-)'"+"'"'ri„<5jv_|i|-2iH,«4 ; (34b) 

A^nAr(^-)'"+""'+'ri„i;A^ = A^(5-)'"+""'+'ri„(5^„|i|_2|m|-i,nfA^ • (34c) 

The passage from the left to the right hand sides has been obtained by moving all the projectors to the right 
with ((25l) . The Kronecker symbols come from nrlln = Sr.rS^n- In (j34ap . the diagonal matrix S~ = 1 — 110 
can be replaced by 1 for n ^ 0. The left hand sides of (IMl) contain the same matrix (5")^ "n„ whose all 
elements vanish but one: ((5"^)^ "n„)jy n ~ ^' W^*^^ acting on a vector a, only a„ gets involved together with 
the last diagonal element Au,jv (Au,)^ ^ of A^,. Using (A^^)^ ^ — N — n and (A-^)^ ^ = n, we eventually 
obtain a vector whose {0, 1, . . . , — 1} components vanish. The last one being 

Tn\,n^x^n ~l~ '^A''— |/ 1 — 2|m| — l,n {N — n)van)- (35) 

words w—{l,m) 
0<|/| + |m|^Af 

The presence of J's forces the relevant words we some on keep |^| + 2|m| ^ A^ — 1 since n ^ 1. Now we have 

Aw,N = a^^N dl'v'^^di'v"'^ ■ ■ ■ (36) 

with the numerical coefficient 

_ ( ,.W+\rni (N - ^i)! (^1 + 2mi - 1)!! {N ^ h ~ 2m, - 1)!! 
""^■^-^ ^> m (?i-l)!! {N-h-iy.\ 

{N-h- 2mi - ^2)! (^1 + 2mi + h + 2to2 - 1)!! (A^ - - 2mi - h - 2m2 - 1)!! 



N-h-2mi\ (Zi+2mi+;2-l)!! (A^ - - 2toi - ^2 - 1)!! 



(37) 



is obtained because (Ai '''')n,n = 1/{N - k + I) and (A[ '" V.w = [k - l)/(iV - fc + 1) for ^ fc > A'. We 
have also introduced the usual notation for any integer n 



n!! = 



2P p! if n = 2p with p integer; 

(2p+iy. ^^^^ 

; — if 71 = 2b + 1 with p integer. 

2Pp\ f & 

I have not been able to prove diectly that the right hand side of pSI) vanishes for an arbitrary A' and for any 
n £ {!,..., A^} but since we are now left with an explicit expression for a linear differential operator acting 
on a„, it can be checked, possibly with the help of symbolic computations, for small specific A^. It is rather 
straightforward to prove that (/i^Z)n„a)^ vanishes for n — N and n = A^ — 1 with arbitrary A^ and that it 
also vanishes for N ^ A and any n e {1, . . . , A"}. 
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1.4 Open questions for possible generalizations remaining with d = 1 



It would be interesting to see if the tools introduced above may be adapted for higher order linear differential 
equations. More challenging, and this will constitute the object future research, we may obtain some interesting 
information on the solutions of some non-linear equations like ([ij . The price to pay is that the non- linear term 
prevents us to keep working in one Vn space only; then, we should work in the whole Vn and probably 
extend condition (i) by factorising a strictly positive polynomial in rather than considering only ip^ . 

2 Higher dimensional cases 

One natural generalization to dimension d is to consider the Schrodinger equation 

d 

A,^ = ^a2^ = 2(V--£;)^ = i;^ (39) 

where x = (x'^)^, G R'^ (we keep definition (flUl) '). In the following, Greek indices always label the dimension and 
we will follow the usual convention of letting implicit the sum from 1 to d over repeated Greek indices imless 

the opposite is specified. The partial derivative with respect to the coordinate is denoted by 9^ =^ d/dx'^. 

We will work within the space VN,d 'Pn of smooth real functions (a„)„ where n = (ni, . . . , Ud) is a 

multi-index with G {0, . . . , N} from which we can construct the set VN,d of functions built from homogeneous 
polynomials, namely having the form 

P4d^{x),vix),x) l5^^^a„(5^)^->^(i-'^)+H (40) 

n 

where \n\ — ni + - ■ - + 71^., {d(p)'^~"' stands for (9i(/7)^~"i (92i/3)^~"^ • • • {ddf)^~"''^ ■ To apply the same reasoning 
that led to gaps in the spectrum, condition (i) will be extended in d dimension by looking for a current J^r = 
{J%)^ whose divergence can be factorised by a (^-independent function times a positive function. More precisely, 
each is associated with an element of VN,d and is constructed in order to fulfill 

D^(j^^{d^{x),^{x),x,E))^^^FN{x,E) . (41) 

Then by integrating it on the whole R**, provided that 

(ii)^ / D^r^d'^x^ f J^d^-'a^^O (42) 

where V is a closed radius whose typical length R tend to infinity (d'^^^tr^ is the measure on its boundary dV 
whose surface growths algebraically with R, therefore any exponential decrease of Jjv will guarantee (ii)(^), the 
condition that Fn{x,E) should change its sign for even TV will hopefully lead to some constraints on E. The 
total derivative is defined as the linear operator in VN,d 

D,[Pa{d^{x),^ix), x)) -^^^^ dl^ip + ^d,ip + d,Pa . (43) 

When d > 1 , we cannot get rid of the second derivatives of ip as easily as for d = 1 because ([5^ provides us with 
only one substitution rulc0: it is only when grouped into a Laplacian, that the substitution X]^=i ^^V' ~ 

^Without further information on V. In the non-generic case of a separable potential, there are in fact d independent substitution 
rules. The substitution rules may also implement some symmetries if there are any, like in the case of integrable systems. 
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be done. If we start looking for a Jn from a generic a, grouping the second derivative in D^J'^ into Laplacians 
will eventually impose some relations on the (a„)„. Following what we have explained in the previous section, 
we will however systematically work up to a total derivative. We can also extend the factorisation in (i)^ to 
other functions f and its derivative whose sign is fixed. Rather than ip^ , we still can apply the argument if we 
manage to obtain 

W'd D^(j';,{d^{x),^{x),x,E)) ^ (B„,^{d^{x),^{x),x,E))\N{x,E) . (44) 
for even N where i?Ar/2 is an homogeneous polynomial in {dip, ip) of degree N/2. 



2.1 Attempt for = 2 

Let us tentatively start with 

J2 = 9^'' ^pdpip + d^Lp + g!^ (45) 

where the g's are smooth functions of (p, x and E. We will start with d functions of type g2, (P functions of 
type gi and (P{d+l)/2 functions of type go such that, without loss of generality, 

gr = 9r- (46) 
All the terms in D^j!^ involving a second derivative in ip can be collected in 

i2grd,p + gr)dl^- (47) 
To construct a Laplacian, we must impose the parenthesis to be anti-symmetric when ^ ^ u : 

gr = ~9r ip + ; (48) 

sr - -97 iP^^), (49) 
and independent of /i when p = ^, that is there are d functions /iq and one function hi such that 

gl^fj-p _ j^P ^-j^Q summation on fi) (50) 

and 

gr^'^ — hi (no summation on /i) . (51) 

Combining (gH) with & leads to .g^"" - g^f"" = -g^^" = -g^J"'' = 50""^ = ffo"^" = -g^"" and therefore 5^"" = 
when {ji,v,p) are pairwise distinct. Collecting in D^j^j!^ the cubic terms in dip leads to 

d 

idv9o'''') d^ipd^ipdpip = (S^/ig) (9p(/3 ^(9^(^)^ . (52) 



Therefore we will take 
Quadratic terms in dip 

d d d d 



d^hP = . (53) 



idp9Q'"' + 9^97) d^(pd^ ^ = ^Y1 ^^^0 d^ipd„ ^ + ^"^0 df^pd^ <^ + X! [('9m<^)^(5¥'^i - X! ^"^"^^ 



fl—1 U — l — 1 V—1 jJ, — l v—1 



(54) 
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Canceling each independent term requires 

d^K + duK = iti^y) (55) 

and 

d 

9fihQ — djyhQ + d^hi = (no summation on /z) . (56) 
i/=i 

The last equation appears as a linear system of d equations that can be rewritten with the help of the d x d 
matrix 2 ~ Ad where {Ad)^,i, — 1. Then, det(2 — Ad) = {—2)'^^^{d — 2) and therefore when d 7^ 2, it can be 
inverted and leads to ^ 

dfih^ = - — - diphi (no summation on /i, d ^ 2) . (57) 

Now Hq does not depend on ip, by (|53p . and then d^hi neither, hence there exist two if independent functions hi 
and hi such that 

hi = hiip + hi . (58) 
The relation (|57p implies (with now an implicit summation on fi) 



hi = ^d^K ■ (59) 



In the special case d = 2, (l56)l leads to 

dihl = d^hl (60) 

and 

d^hi = . (61) 



Then we can keep (j58p together with (j59l) even for d = 2. 
Collecting all the previous relations we get 

D^J!^ = {2h'^^d ^ + d, g"/ + d^gt^)d^^ + /iiA^ ip + 9^52 (62) 

where of course, we can use the substitution (15^ . Without loss of generality we can take g'^^ = for /i 7^ by 
possibly redefining 

g^iif,x,E)^g^ip,x,E)-y2 g\^{^' ,x,EW (63) 

since d^gi^ ^ d^gi^ — X]'5=i Jn d^d^gi^' {f' , x, E)d(p' = d^^g!^ because the integrand cancels by P5|) . To cancel 
the parenthesis in (j62p we must take 

9^92 = -9f_,hi - 2h'^v ip (64) 



The dependence in appears only through ((58)) and we can immediately integrate the last relation 

9^2 = - (^5^ hi + /i^t') - (9^ ~hi)V + 32^ (65) 



where 32 is (^9- independent. With this expression, (p^ becomes 



-Dm>^2 = -V' 



-Arfft-i - vhi + 9^(wft.o) 



- (y5[Ad/ll - w/ll] + 9^x32 ■ (66) 
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The last term is irrelevant because it is a total divergence and can be reabsorbed in the definition of J2 . The 

second term does not contribute also since by integration by part it can be converted to hi[/S.d ~ v\lp which 
vanishes, with the use of (|59p. the first term can be further simplified in order to keep /iq only. 
To sum up, condition (i)^ can be obtained for N ^ 2 with 

F2 - ^ Add,,h^ ~^vd,X~ h^df^v (67) 

with /ig being any d smooth functions such that 

dihl = d2hl^--- = ddK; (68a) 
d^^K - -duK (m^'^)- (68b) 

Our freedom of choosing J2 has therefore being reduced first because eliminating the second derivatives of Lp 
through its Laplacian impose severe constraints and second because, as in the d = 1 case, many different initial 
choices lead to the same F2] in others words the linear application from a to has a non-zero kernel. 

For d = 1, we have one test function ho — go ~ oq at our disposal but without any restriction on its 
derivative and we can easily checked that, for d = 1, ([57)) is (|3b[) . 

For d = 2, we have F2 = —d^{vhQ) and surprisingly (j68l) appears to be the Cauchy-Riemann conditions 
for + i/iQ to be analytic. Anyway, if we compute 

F2dxdy^ [ v^h-da (69) 
s Jds 

on a surface S G whose boundary dS coincide with the energy level V{x, y) — E, that is u = 0, the integrand 
of the right hand side vanishes (da is an infinitesimal 2(i- vector normal to the curve dS pointing outwards say) . 
Therefore whatever choice we make for ho, for any energy E belonging to the image of V (where we know the 
spectrum lies), F2 changes its sign and no information can be obtained further. 

For d ^ 3, the constraints (|68p are so strong that they limit the choice of ft-o's to polynomials in x of 
degree at most 3. Indeed all the third derivatives of ho must cancel (up to equation (|73p included, any pair of 
distinct Greek letters denote any pair of distinct values of indices and no summation over repeated indices is 
left impliciiEl). First, 

dl^K = -dLK = dl,K = -dl,K (70) 

and therefore df^^^h^ — 0; then d^^^ph^ ~ 
Furthermore, 

dl,.,.K = -dlpX - = (72) 

since for d > 3 we can always find an index p distinct from both fi and i^; eventually we have 

'^fj.fj.fj."'Q — '^fj.fj.u"'Q — '^p,uu"-0 ~ '^pui^"-0 ~ '^ppu"'0 ~ '^ppp"'a ~ '-' pti.ti."'0 ~ '^ppp"'0 V ' "^z* 

and all these third derivatives actually vanish as well. Now /io(x) being a polynomial of degree at most two 
in X, the constraints (1551) on its coefficients leads to the general form 

h^{x) = h^{0) + kx" + A^^x"" - i l^x^ +x''l-x (74) 



^1 am grateful to Olcg Lisovvvi4l for providing the following arguments that concisely and rigourously proved my first guess of 
(El- 
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where A is a d x d constant antisymmetric real matrix, k a real constant, I a real constant d- vector ; the 
Cartesian product I ■ x — If^x^ is used. Then, from (|67p we get: 

F^ix, E, A) = 4(S - V{x)) {k + l-x)~ (/i^(0) + kx^' + A^'^x" ~ ^ Fx^ + x^" I ■ x)d^j,V{x) . (75) 

Unlike what occurs for d = 1 where we are free to construct F2 from a whole set of test functions x i~> 00(2^), 
for d > 3 we are left with only [d^ + Sd + 2)/2 free x-independent parameters, namely A — (/iq (0), k, A^^, Z'^) . 

Now the boundaries of the gaps must belong to the solutions of The linearity of F2 in A simplify 
considerably the computations. The conditions (pc)) imply (I5al) and are equivalent to 

E = V{x) ■ = . (76) 

Condition (I6b[) is guaranteed if we choose for instance 

Ki^^^kx^^A^^.x^-,\l^xl^x^l. ill) 

and condition (pd) is generically fulfilled. 

Therefore, with our method, possible candidates for the gap boundaries are the critical points Xc of V 
which is not a surprise from a semi-classical point of view. With this method we cannot expect to find more 
interesting and more relevant piece of information. Actually, some inequalities concerning the global spectrum 
may be obtained if we are to maintain the sign of (j75p . specially once a specific V is given; but our ambition 
was, as we have shown in [2] for d = 1, to obtain some local information in the very core of the spectrum. 



3 Ending remarks 

3.1 Simplified starting point 

For N = 2 and any d we have shown directly that without loss of generality we could have started with no term 
in ip in J2 , that is with a current such that 

d^J^{d^{x),'f{x),x,E) ^0 . (78) 

Indeed, the are ip independent, see (j53|) . and we could have taken hi = with no consequence on the 
result (pT)) . For d = 1 and any N, we have proven this result through p^ : we obtain all the possible F^v's even 
if we restrict our self to a = (oq, 0, • • • ). This result can be also obtained for any d and any N in another way. 
First remark that if we start with J'^(^dip{x),(p{x),x,E'j the second derivatives in 

Z?^ Jj^ = Jj^ + d^J'^ df,Lp + da^^J'jij df,„ip (79) 

can be eliminated with the help of ([55]) if and only if 

da^^J^ = LS^i + T^^. (80) 

where S is the Kronecker symbol, W(^d(p{x), ip{x), x, E) an antisymmetric d x d matrix and L(d(p{x), (p{x), x, E) 
a function. Then after the substitution of p9p . an integration by part can be made and we have 

D^^J^ = d^J^ + Lvp - D^id^J^)^ + D^id^J^) . (81) 

By computing D^j^{dipJ'^) in the same way, and iterating the process up to infinity, we find that D^J'^ can be 
written like 

D,r^ = v^Y. + E ^-^d:^d,r^ + d, (r^^ f; ^-^d;A ; . (82) 

ri=0 ■ n=0 ' V n=0 ' / 



13 



Up to a total derivative, starting from any J^r, we therefore are always led to 



where Hq in the evaluation at (p — that can be expressed as 



dJUoJ^) = vipUoL + Uod^J^ (83) 



no-E^^^." (84) 



n 

n=0 



for any analytic function of (p. Therefore even if we start with a Jn whose component takes the general form 
working up to divergence terms, we will be led to the same identity as if we had started with all the a„ such 
that iV(l — 0?) + |n| > being zero (this is consistent with (|T^ ). In = 2, all the g's in (|^5|) could have been 
taken independent of ip from the very beginning. 

3.2 Remark about condition (i)'^ 
For d ^ 2 and N = 2, condition (i)^ with 

Bi = bfid^ip + boip (85) 

where (bi = (&i )^, &o) are d+1 smooth real functions of x, does not provide F2 explicitly but rather leads to a 
linear differential equation for it. For instance, when d = 1 we get 

i blF^' + i2b,b[ - b,b,)F.^ + [ib,b[Y - vbl - {b.boY + bl]F2 = ia^" - 2va'o - v'ao (86) 

which reduces to (j3bp when we make the simplest choice bo — 1 and 61 = 0. For d > 1, by reproducing the 
same line of reasoning as in § 12.11 we get a second order linear partial differential equation for F2 



{b^YAdF^ + b ^1 • 9^.h - b^bo d^F2 



Adibif ~ -vibif - d^ibf^bo) + bl 



Fo 



2 — d 2 

Add^h'^ - -V d^h'^ - hl^d^,v . (87) 



instead of (p7| obtained for 60 = 1 and 61 = 0. The constraints on /ig now involve the 6's and are entangled 
with F2: 

dihl~libl)'F2 = d2hl-lib',fF2 = ---^ddh^o'libfrF2; (88a) 
+ d,h^, = b'^b'(F2 ilJi^v). (88b) 

For d > 1, the only way to get rid of F2 from (|88ap is to take b^ independent of /i, and eventually 61 = if we 
want (|88bp not to involve F2 either. Then, since F2 and 69^2 have the same sign, ((87|) take us back to ([67)) that 
is to case (i)^ where 60 = 1. 

Even for d = 1 and a specific w, I did not exploit further these possibilities, but the choice of 6's for which 
can be solved explicitly is rather limited (not to speak of the control of the sign of its solutions). In any case of 
course, we ought to work with simpler differential equations than the Schrodinger equation itself ! 

3.3 Attempt for d = 2,N = A 

The second remark concerns an attempt to obtain (i)^ for d = 2 and iV = 4. Using the same ideas as in the 
case N — 2 and with 

4 

Jit = E 54-™.™ (^1^)'"" (^2^)" (89) 

n=0 
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with, according to the last remark, taking g as (p- independent. The constraints (|80p imposed by the ehmination 
of the second derivatives of in D^J^l lead to two independent functions instead of the ten g's. But if we go 
further to eliminate the cubic terms in dip, these functions must vanish identically and no non zero g can be 
found this way. Extending (i)^ to (i)^ with 

B2 = b2o {diip)"^ + hii diifi 82 (p + bo2 {d2^Y +biQipdiip + 601 ipd2^p + b^Qip^ (90) 

where the &'s are smooth function of x, leads to the same conclusion. 
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